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Abstract. During NASA’s GLObal Backscatter Experiment (GLOBE) II flight mission over the 
Pacific Ocean in May-June 1990, extensive aerosol backscatter data sets from two continuous 
wave, focused CO 2 Doppler lidars and an aerosol microphysics data set from a laser optical particle 
counter (LOPC) were obtained. Changes in aerosol loading in various air masses with associated 
changes in chemical composition, from sulfuric acid and sulfates to dustlike crustal material, sig- 
nificantly affected aerosol backscatter, causing variation of about 3 to 4 orders of magnitude. Some 
of the significant backscatter features encountered in different air masses were the low backscatter 
in subtropical air with even lower values in the tropics near the Intertropical Convergence Zone 
(ITCZ), highly variable backscatter in the ITCZ, mid-tropospheric aerosol backscatter background 
mode, and high backscatter in an Asian dust plume off the Japanese coast. Differences in aerosol 
composition and backscatter for northern and southern hemisphere also were observed. Using the 
LOPC measurements of physical and chemical aerosol properties, we determined the complex 
refractive index from three different aerosol mixture models to calculate backscatter. These values 
provided a well-defined envelope of modeled backscatter for various atmospheric conditions, giving 
good agreement with the lidar data over a horizontal sampling of -18,000 km in the mid- 
troposphere. 


1. Introduction 

Aerosols form a complex system with varying composi- 
tion, size, number distribution, morphology, and shape. 
These parameters govern the aerosol scattering properties at a 
given electromagnetic radiation wavelength and in turn are 
affected by the meteorological conditions, altitude, and 
geographic location of the air mass. As an example of one of 
these parameters, aerosols can occur in various compositions 
[Prospero et ai, 1983; Clarke and Porter , 1991; Clarke, 1993] 
such as sulfuric acid, sulfates, nitrates, minerals, dustlike crus- 
tal material, and sea salt, which can exist in pure form, as 
composites, or in complex mixed phases [Tang and 
Munkelwitz , 1977; Tang et ai, 1978]. This variety leads to a 
wide range of complex refractive indices of aerosols, thereby 
affecting their scattering. 

It is difficult to measure all aerosol parameters required to 
exactly predict the scattering, since most of them cannot be 
simultaneously determined by one instrument. A combination 
of instruments is generally required to determine various 
aspects of aerosol microphysics [Willeke and Baron , 1993]. 
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In addition, some parameters that are very difficult to measure, 
such as shape or morphology, are usually approximated 
through the use of different aerosol models. However, these 
models may lead to different scattering results [Srivastava, 
1988]. To ascertain a model's validity, scope, and limitations, 
model results need to be compared with direct scattering meas- 
urements. Further, such comparisons of modeled scattering 
using in situ aerosol microphysics data inputs can lead to 
better understanding of aerosol chemical properties and mor- 
phology of the complex aerosol systems and how they can 
affect their scattering. 

With increasing use of lidar remote sensing of various 
atmospheric parameters, more detailed aerosol backscatter 
variability is being detected [Post, 1984; Trait and Menzies , 
1994; Vaughan et ai, 1995; Rothermel et ai, 1989, 1996a], 
enhancing the need to better model aerosol backscatter for 
various applications, e.g., winds [Baker et ai, 1995]. 
Airborne and ground-based, pulsed and continuous wave (CW) 
lidars provide large-scale remotely sensed aerosol backscatter 
data [Rothermel et ai, 1996a; Menzies and Tratt, 1994]. CW 
CO 2 lidars provide data with very high sensitivity and high 
spatial/temporal resolution even in clean conditions 
[Srivastava et ai, 1995], measuring backscatter as low as 3 x 
10' 11 m" 1 sr" 1 in aerosol loading of fractions of a microgram 
per cubic meter of air. On the other hand, advanced aerosol 
counters can provide both detailed aerosol size distribution 
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and compositional information on relatively short timescales 
[ Clarke , 1991, 1993]. With the use of the in situ aerosol 
microphysics data, aerosol backscatter from different compo- 
sitions can be appropriately modeled. These modeled aerosol 
backscatter results can be compared with direct, coregistered 
lidar backscatter measurements, providing verification of the 
modeling of atmospheric aerosol scattering. As the ambient 
atmospheric conditions change, aerosols can undergo chemi- 
cal and physical changes, giving rise to complex mixtures, 
possibly requiring different models. Thus airborne measure- 
ments with aerosol counters and lidars provide detailed varia- 
tions in aerosol microphysics and scattering, respectively, 
while their intercomparison can indicate variations in aerosol 
mixtures and the best models to use. 

2. Measurements 

To extend the understanding of atmospheric aerosol back- 
scatter distribution, the National Aeronautics and Space 
Administration (NASA) sponsored a GLObal Backscatter 
Experiment (GLOBE) flight survey mission over the Pacific 
Ocean in 1989 and 1990 [Bowdle et ai, 1991]. GLOBE con- 
sisted of a complement of various instruments on NASA’s 
DC-8 aircraft, measuring various atmospheric aerosol proper- 
ties. Among the various instruments flown were two NASA 
CO 2 Doppler lidars measuring aerosol backscatter [Rothermel 
et ai, 1996b]. In addition, one of the aerosol counters flown 
was the University of Hawaii’s preconditioned laser optical 
particle counter (LOPC) [ Clarke , 1991, 1993], providing in 
situ measurements of microphysical and chemical properties 
of the aerosols. The measurements were identified by the 
universal time coordinated (UTC) code at the time of data 
storage. 

2.1. Lidar Backscatter Data 

Direct lidar backscatter measurements were obtained with 
NASA/MSFC’s focused, CW, coherent, homodyne CO 2 Dop- 
pler lidars operating at 9.1 and 10.6 pm wavelengths. These 
lidars transmitted focused radiation at -50 m away from the air- 
craft in front of the aircraft’s left wing into the atmosphere and 
received Doppler-shifted backscattered signal from the aero- 
sols in the lidar sample volume. The measured signal spectrum 
was converted to signal-to-noise ratio (SNR) data [Rothermel 
et ai , 1996b; Vaughan et ai, 1989]. Measured SNR was 
converted to aerosol backscatter, p (m" 1 sr" 1 ), by using 

radiometric lidar calibrations [Jarzemhski et ai , 1996]. 

2 

Coherent detection and a large lidar sample volume of -6x10 
m 3 [Cutten et ai, 1996] allowed high sensitivity and detection 
of low signals for clean aerosol conditions, making it 
possible to acquire good quality high-resolution data with 5 s 
integration time (an aircraft speed of -200 m s" 1 implies 
spatial resolution of -1 km). A larger lidar sample volume 
implies more particles sampled in the lidar beam, leading to 
smoother signal peak averaged over the backscattered return 
from all particles and allows a better estimation of SNR in low 
signal conditions. In clean mid-tropospheric conditions the 
9.1 and 10.6 lidars detected aerosol p as low as 5 x 10‘ 12 m" 1 
sr’ 1 and 2 x 10' 11 m’ 1 sr' 1 , respectively. However, near the 
detection threshold there is more variability in both lidar data 
sets due to the added uncertainty of extracting low signals out 
of noise. The uncertainty for moderate to high P is -24%; 


however, for low p < 5 x 10’ 1U m" 1 sr’ 1 it is -36%. Details of 
lidar instrumentation, calibration, and sensitivity study are 
given elsewhere [Rothermel et ai, 1996b; Jarzembski et ai , 
1996]. 

2.2. LOPC Aerosol Microphysics Data 

The LOPC system [ Clarke , 1991] consists of a modified 
laser aerosol spectrometer (LAS-X, PMS, Boulder, Colorado) 
operating at 0.6328 pm wavelength with 256 size bins for 
particle diameters ranging from 0.15 to 7.0 pm. It has five 
selectable inlets heated to different temperatures such that 
some compositions, depending on their volatility, can evapo- 
rate while the residual composition’s size distribution can be 
measured. The LOPC system has the unique capability of pro- 
viding information on separate size distributions for some of 
the most important aerosol components: sulfuric acid, ammo- 
nium sulfate and bisulfates, and dustlike crustal material. These 
compositions are the ones most frequently encountered in mid- 
tropospheric aerosols, and their fractions can be inferred from 
the LOPC [Clarke and Porter , 1991; Clarke , 1993]. The LOPC 
was mounted on the right side of the DC-8 aircraft in front of 
the right wing. Detailed descriptions of the instrumentation 
and calibration are given elsewhere [ Clarke , 1991; Porter et 
ai, 1992]. 

The LOPC’s preheater system successively heats the aerosol 
sample, drawn in from the ambient atmospheric conditions, to 
40°C, 1 50°C, and 320°C in channels Chi, Ch2, and Ch3, 
respectively. Aerosols are heated to 40°C in Chi, giving a 
nearly dry total aerosol size distribution at relative humidity 
(RH) of -10% within the LOPC due to partial evaporation of 
water component. In Ch2, aerosols are heated to 150°C, where 
the volatile sulfuric acid (SA) component of aerosol evapo- 
rates. In Ch3 they are further heated to 320°C, where the 
ammonium sulfate (AMS) component and any bisulfates 
evaporate, leaving the residual nonvolatile, refractory, dust- 
like (DS) crustal component. Differentiation between AMS and 
bisulfates is not possible. Each channel gives the number of 
particles in the given size bin as dN/dD, which can be inte- 
grated to give the total number of particles No [particles/cm 3 
of air] in that channel, indicative of the aerosol loading. To 
improve count statistics for low concentration of coarse aero- 
sols and shorten accumulation times, a virtual impactor (VI) 
[ Clarke , 1991] is used along with the LOPC in channels Ch4 
and Ch5. The VI Ch4-measured aerosol size distribution heated 
only to 40°C gives total dry aerosol content, while VI Ch5 
heated to 320°C gives the distribution of the residual 
refractory dustlike crustal component after both the SA and the 
AMS components have evaporated. 

In general, the uncertainty of the LOPC data is roughly 10% 
to 20%. However, for clean conditions with fewer large- 
particle counts the uncertainty can be higher depending on the 
sample time. In these cases the sample time was usually made 
longer to minimize this uncertainty. LOPC size distributions 
and forward scattering spectrometer probe (FSSP) data have 
been compared, giving fairly good agreement [Porter 
et ai, 1992; Cutten et ai, 1996]. 

The information available from the LOPC channels was used 
to model aerosol /?. However, for comparison of the lidar p and 
modeling using the LOPC data it is important to note certain 
key differences in the two data sets. LOPC data averaging, as it 
cycles through its five channels, can vary from approximately 
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1 to 3 min per channel, leading to about 5 to 15 min per 
observation, implying spatial resolution of ~60 to 180 km. 
During an LOPC observation there can be nearly 60 to 180 
lidar data points exhibiting the fine-scale aerosol variation 
within the LOPC sample. Air mass changes occurring during 
an LOPC observation are averaged over the finer-scale varia- 
tion shown by the higher-resolution lidar data. In cloudy air or 
near cloud boundaries the aerosol loading can change so rap- 
idly such that modeling with the LOPC data (averaged over 
highly variable conditions) can have limitations due to both 
long sampling times and loss of water in heated channels, 
which cannot be easily accounted for at high RH (discussed in 
section 3.2). Furthermore, the LOPC sample volume of ~1.7 x 
10” ? m 3 is several orders of magnitude smaller than the lidar 
sample volume of -600 nr . In clean conditions, where there 
are only few large aerosols, even with the inclusion of VI 
channels LOPC can underestimate their number concentration 
and the modeled aerosol optical properties. Lidar, with its 
higher sample volume, is less susceptible to lower large parti- 
cle count statistics and can perform well as long as the signal 
is above its detection threshold. 


3. LOPC Data Analyses 

The aerosol microphysics data provided by the LOPC must 
be interpreted for the aerosol /J modeling to be performed. The 
volume fraction, v, of each component in each size bin was 
estimated from the distributions measured in different chan- 
nels. In addition, LOPC operates at aircraft cabin temperature, 
resulting in dry aerosol size distribution well above the ambi- 
ent. Hence the aerosol diameter, D, provided by LOPC was cor- 
rected to ambient conditions to get the actual aerosol diameter, 
D^ y compensating for possible loss of water content in the 
heating process (discussed in section 3.2). 


3.1. Volume Fractions of Various Aerosol Components 

By using the information in each channel the various com- 
ponents can be differentiated, and their fraction can be esti- 
mated. Thus, for a given component the differential volume 
fraction, v, (where j represents either SA, AMS, or DS), in 
each size bin can be determined from the difference of the 
distributions in the various channels. For channels Chi, Ch2, 
and Ch3 these are given as 
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and for the VI channels they are given as 
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where dN/dD represents number of particles per cubic centime- 
ter of air in a given size bin. Each component’s v r - can be 
either integrated over all bins to get total volume fraction, 
V? , or converted to differential volume in each bin and then 
integrated to get the component’s total volume, V, . 

Figure 1 shows some aerosol number distributions, dN as a 
function of dry particle diameter D, measured by LOPC in 
channels Chl-Ch5. Channels Chl-Ch3 are shown in Figures 
la and lc, while VI channels Ch4 and Ch5 are shown in Fig- 
ures lb and Id. The distribution measured in Chi (solid line) 
represents nearly dry total aerosol content at 40°C. When the 
aerosol is heated to 150°C in Ch2, SA content and any residual 
water volatilizes, giving a different number distribution 
(dashed line). In Ch3 the aerosol is further heated to 320°C, 
where the sulfate component evaporates, leaving a dustlike 
refractory component distribution (dotted line). The difference 
between the distribution at 40°C and at 150°C gives an esti- 
mate of the concentration of SA content in the aerosol 
(equation (1)). Similarly, the difference between distribution at 
150°C and 320°C provides the AMS content of the aerosol 
(equation (2)). The residual nonvolatile DS-like crustal or 
refractory content of aerosol is given by the distribution 
obtained at 320°C (equation (3)). Inclusion of VI in channels 
Ch4 and Ch5 shows enhanced retrieval of large particle 
counts. Since the VI channels do not include a middle channel 
heated to 150°C, the difference between Ch4 and Ch5 gives the 
combined fraction of the SA and AMS component (equation 
(4)), while Ch5 gives again the DS component (equation (5)). 
Thus, to distinguish between SA and AMS from channels Ch4 
and Ch5, the ratio of the fraction of SA to AMS from the 
difference of channels Chi, Ch2, and Ch3 (equations (1), (2), 
and (3)) is used. This allows the use of more robust Vl-included 
size distribution measurement in channels Ch4 and Ch5 with 
component information available from channels Chl-Ch3. 

The example of LOPC size distribution data shown in Figure 
1 is from GLOBE II, flight 12, in transit from Darwin, 
Australia, to Tokyo, Japan, obtained at -8 km altitude. 
(GLOBE flight track details are given by Cutten et al. [1996]). 
Figures la and lb show low aerosol concentration in the rela- 
tively clean subtropical conditions encountered in the early 
part of the flight around - 1 0°N. Channels Chl-Ch3 show 
mostly SA and some AMS component with little or no DS. 
Mostly SA and AMS were found in subtropical conditions with 
some fine DS fraction as possible residue of continental effect. 
Similarly, Figures lc and Id show size distributions in Chl- 
Ch5 with higher aerosol loading due to an Asian dust plume 
event mixed with some anthropogenic aerosols near the 
Japanese coast in the last part of the flight. The difference 
between the channels again gives an estimate of the volume 
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Figure 1 . Aerosol size distribution data measured by LOPC in five heated channels, Chl-Ch5. Two different 
data sets from GLOBE II flight FI 2, on May 31, 1990, from Darwin, Australia, to Tokyo, Japan, at ~8 km 
altitude show examples of: (a, b) clean tropical/subtropical aerosol loading at about (7°N, 134°W), and (c, d) 
heavier aerosol loading at about (34°N, 139°W) in an Asian dust plume. 


fraction of various components of the sampled aerosols. In 
this case, larger particles were encountered with a significant 
amount of DS and some AMS and SA component. These figures 
demonstrate the changes in size distributions due to volatility 
of the aerosol components in different channels at given tem- 
peratures, as well as changes in the aerosol components under 
different atmospheric conditions. 

3.2. Correction to Ambient Aerosol Size 

Because of the large temperature increase in bringing free 
tropospheric ambient aerosols into the aircraft the RH drops 
significantly, causing the aerosols to lose most of their water 
content. Therefore the size obtained is nearly a dry diameter, 
P, measurement at the LOPC channel temperature and RH 
rather than at the ambient temperature and RH. The actual 
aerosol diameter, D A , has to be determined at ambient atmos- 
pheric conditions outside the aircraft. Hence measured D was 
corrected for ambient aerosol water content to get D A , which 
was used for modeling /?, The procedure for extrapolating to D A 
was done in postprocessing of the data by determining the 
growth of the aerosols from LOPC RH back to ambient RH, 
both of which were recorded as part of the data set. A growth 
factor, G, was applied to the measured distribution to get the 
consequent aerosol size change, resulting in a near-ambient 
aerosol size distribution. 

Different aerosol compositions grow differently on the 
basis of their hygroscopicity. For example, SA, being very 


hygroscopic, changes size dramatically as RH is changed, 
whereas DS-like crustal material hardly grows at low RH (for 
example, in the middle and upper troposphere). Growth for 
each component, i, of the aerosol can be given in terms of a 
growth factor Gi=(D A /D) it where D and D A are the dry and ambi- 
ent diameter, respectively. Growth curves given by Tang and 
Munkelwitz [1977] and Tang et al. [1978] for SA and AMS 
were used for estimating growth factors G SA and G AMS . How- 
ever, growth for DS, G DS , can be very complex [Hand, 1976; 
Hanel and Lehmann , 1981], depending on its mineral content 
and possible sulfate coatings; hence for simplicity an average 
activity factor of -0.1 has been assumed to get Gq$. Growth is 
primarily dominated by the SA component, as it is signifi- 
cantly more hygroscopic than AMS or DS. Once the volume 
fractions of the components have been estimated by using 
equations (l)-(5), G, can be determined for each bin size. 

The method of estimating aerosol growth also depends on 
the type of aerosol mixture. In the simple case in which each 
aerosol is made of a pure single component with no mixing, 
D a for the individual aerosol component for the each size bin 
can be estimated by multiplying D ( by its respective G 4 . In 
cases in which aerosols can have two or more mixed compo- 
nents, growth would depend on their v ( in the mixture; there- 
fore a volume-weighted approach was used. The effective 
volume-weighted growth, G e , for mixed aerosols can be 
obtained from the component’s G, on the basis of its v, in each 
size bin as 
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G e ={l(G,.) 3 *v ( } l/3 , (6) 

where the summation over i stands for the components present 
in mixed aerosol. In the case of =1, (6) reduces to G e = G,, 
going back to the single-component growth. The ambient 
mixed aerosol diameter is given by D A = G e * D. The G e 
changes slightly for each size bin, as the composition frac- 
tions can change in different size bins, as seen in Figure 1. For 
a complex aerosol system in which composition can vary for 
different sizes this size dependent growth provides a simple 
but more realistic estimation of D A in comparison with assum- 
ing a single growth factor for the entire distribution. In addi- 
tion, it allows for two or more components to be mixed, 
depending on different aerosol models. 

The growth factor was further used to estimate the volume 
fraction of the water content lost on heating the aerosols in 
the LOPC channels in each size bin. The volume fraction of 
water for each component, v Wti , can be given in terms of the 
growth factor G { for the pure aerosol case as 

1 


where for the mixed aerosol case the subscript i is replaced by 
e , using G e for mixed components. The water volume fraction 
contribution to backscatter can be taken into account by using 
v w estimates. 

Even though G SA dominated the aerosol growth, a volume- 
weighted growth approximation was adopted for mixed aero- 
sols to account for any AMS and DS effect and for consistency 
as a similar approach was employed in the commonly used 
volume-weighted effective refractive index discussed next. 
Both the Gj and the G e may be limited for very high RH 
(beyond the deliquescent point of solids) but are fairly adequate 
for size changes for p applications in the range of RH of 5% to 
60% in the middle and upper troposphere encountered during 
GLOBE. 


4 • Aerosol Backscatter Modeling 

Atmospheric aerosol backscatter, p, at a given radiation 
wavelength, A, depends on D A and composition characterized 
by the complex refractive index, m. For different v,-, m can be 
estimated from LOPC data (equations (l)-(5)). Various aerosol 
mixing models can be used, since atmospheric aerosols are 
continuously changing and interacting as the ambient condi- 
tions are changing, giving rise to complex structures. Conse- 
quently, aerosol compositions can occur in a variety of con- 
figurations or mixtures. It is almost impossible to predict 
exactly how these compositions are mixed in an atmospheric 
aerosol ensemble; however, some plausible models can be 
formulated for certain conditions to simulate them. In very 
clean conditions frequently encountered in remote marine 
atmosphere, newly produced aerosols may exist separately as 
pure components for short periods of time (a few minutes to a 
few hours). As these aerosols age and mix in changing air 
masses, they tend to lose their pure and simple configuration, 
becoming partially mixed because of growth by diffusion, 
deposition, coagulation, or cloud processes. They may be 
more stable and have longer lifetimes (of the order several 
hours to days). In marine mid-troposphere these may be pre- 
dominantly acidic sulfates with possibly some fine dust com- 


ponent residual of continental effect, or on a few occasions 
they may contain some large sea-salt particles in very strong 
convective updrafts. The dust component may or may not 
become composited with sulfuric acid and sulfates. In aged 
continental air masses there can be significant increase in dust 
or pollution-type aerosols. These can have varying lifetimes 
(a few hours to days or even weeks) and can be advected over 
considerably long distances, thus giving rise to a high 
probability of completely mixing with acids and sulfates in 
internally mixed composites or even coated particles. 

To deal with these different scenarios, three simple aerosol 
models were chosen to calculate m and p to assess the effect of 
different mixing: (a) pure, (b) partially mixed, and (c) com- 
pletely mixed components. These three models provide a 
fairly good range of possibilities for various intermediate 
cases that may exist. Figure 2 shows a simplified schematic of 
the three models with random volumes of different 
components. In the mixed cases the volume fractions are 
represented as pie fractions, which would, however, be 
randomly dispersed within the particle [ Srivastava , 1988]. 
Another specific model deals with an intermediate case of 
coated particles that may exist in the atmosphere. This model 
is more difficult to evaluate, since it is hard to accurately 
determine the coating thickness or its chemical properties if 
the coating is a complex mixture. Thus it would require 
additional assumptions on the size and kind of thickness; 


I 1 SA 




DS 


(a) External Mixture (EM) 




(b) Mixed Phase Sulfate (MPS) 







(c) Internally Mixed Composite (IMC) 



Figure 2. Schematic of three different aerosol models: (a) 
external mixture (EM), in which all components, SA, AMS, 
and DS are separate, (b) mixed-phase sulfate (MPS), in which 
only SA and AMS are randomly mixed together with varying 
volume fractions while DS is separate, and (c) internally mixed 
composite (IMC), in which SA, AMS, and DS are randomly 
mixed together within each particle in varying volume frac- 
tions. For simplicity in mixed particles the volumes are sche- 
matically shown as pie fractions, which, however, would be 
randomly distributed within the particle. 
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hence for simplicity, only a general internally mixed 
composite model is treated here in case (c). 

Aerosol backscatter modeling was performed by using Mie 
theory [Mie, 1908], assuming a spherical shape for the aero- 
sols. This fits well for SA and also can be applied satisfacto- 
rily for AMS, which, being hygroscopic, may generally be 
slightly wet amorphous structures resembling rough spheres. 
The DS component, however, would be nonspherical, and for 
these cases, Mie theory is an approximation. However, for 
long wavelengths this assumption has given reasonable 
results [Cutten et al., 1996]. For shorter wavelengths, non- 
sphericity may become significant. Complex refractive 
indices of pure SA, AMS, and DS aerosols [Srivastava et al., 
1992, 1995; Kent et al., 1983] and water at the two CO 2 lidar 
wavelengths were obtained from refractive index 
measurements [Shettle and Fenn , 1979; Toon et al., 1976; 
Palmer and Williams , 1975; Hale and Querry , 1973]. 

4.1. External Mixture Model 

In the simple external mixture (EM) aerosol model, SA, 
AMS, and DS are assumed to be separate pure components 
(Figure 2a). This situation may be expected to occur in clean 
mid-troposphere, especially in the tropical regions, where 
new particle formation is more likely [ Clarke , 1993], giving 
less chance for intermixing. In this model therefore the 
effective complex refractive index, m l+ vv* of each aerosol is 
calculated from its pure component fraction along with water 
fraction as estimated by the growth (equations (6) and (7)). 
Thus m r+ vv ambient aerosols can be given by the volume- 
weighted average of each m, with v, and m w with v W i as 

m SA + W - v SA m SA +v W,SA m W (8a) 

m AMS+W = v AMS m AMS+ V W,AMS m W ( 8b ) 

m DS+W “ v DS m DS + v W,DS m W • (8c) 

The water content associated with the AMS and DS component 
was almost negligible for the low mid-tropospheric RH. The 
most significant effect of water is felt by the hygroscopic SA 
component; however, the growth in the radii is slightly offset 
by the dilution effect of water to m SA+w . With m^ w estimation 
at a given wavelength (equations (8a-8c)), the single-particle 
backscatter cross section, <7 j( D A , m i+w , X) (m 2 
sr* 1 ) is calculated using Mie theory. For the EM model the total 
backscatter coefficient, P is given by 



where [dN/dD A ]$ A or AM $ or DS is obtained from channels Chl- 
Ch5 given by equations (l)-(5). 

4.2. Mixed-Phase Sulfate Model 

In the mixed-phase sulfate (MPS) model, SA and AMS are 
mixed internally in one particle and DS component remains 


separate (Figure 2b). These may occur in slightly older, par- 
tially mixed air masses where the probability of mixed-phase 
ammoniated sulfates is high though the DS component may 
not have had a chance to get mixed in. These can chemically 
react and form bisulfates or other intermediate sulfates. How- 
ever, the optical scattering properties of bisulfates are not 
well known in the infrared. Hence a volume-weighted physical 
mixture of the individual SA and AMS properties has been 
assumed as internally mixed, while DS is treated separately as 
an external component. In this case the effective complex 
refractive index of mixed sulfate, m MPS , is calculated for each 
size bin by using complex refractive indices of SA, AMS, and 
water along with their v, as 

m MPS - V SA m SA + V AMS m AMS + v W, MPS m W > ( 1 ^ a ) 

and for DS as given before by 

m DS+W - V DS m DS + V W.DS m W* (10b) 

where the vw.mps iS the water volume fraction in the MPS com- 
ponent (equation (7), with G e = G^pj) and v w DS is in the 
unmixed DS component. Volume fractions of the components 
are obtained from the difference of the LOPC Chi to Ch3 as 
given in equations (l)-(3). As we noted above, 
°MPS K , m ww , A) for the mixed sulfate is calculated with 
the mixed sulfate effective complex refractive index and for the 
DS component separately. For the MPS model the total back- 
scatter, p, is obtained by 



where [dN/dD A ]s A+A \is I s obtained from equation (4) as the 
difference between Ch4 and Ch5. 

4.3. Internally Mixed Composite Model 

This internally mixed composite (IMC) model contains SA, 
AMS, and DS mixed in one composite aerosol (Figure 2c). 
This type of aerosol may occur in fairly aged complex air 
masses, such as pollution or continental dust plumes, that may 
be advected for long distances before dissipating like the 
Asiatic dust plume over the Pacific Ocean encountered off the 
coast of Japan. This would also be suited for complex 
boundary layer aerosols, high-humidity regions, and cloud- 
pumped aerosols. The effective complex refractive index, 
m /wc» °f the composite aerosol is obtained by volume 
averaging of m^and myy as 

m IMC = V SA m SA + V AM m AM + V DS m DS + V W,IMC m W » ( 1 2 ) 

where v WIM ^is water volume fraction in the composite 
aerosol (equation (7), with G e = G/^c). In the IMC model, 
®imc( ®A> m /MC> f° r the composite aerosol is calculated with 
m iMC t0 S* ve total p as 

A , / JdN~ 

P-^IMc{ D A fm iMC^) a » (13) 

- A JSA+AMS+DS 

where [dN/dD A ] SA+AMS+DS is for total aerosol distribution 
obtained from Ch4. 
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5. Comparison of (5 Modeling With Lidar Data 

Results of backscatter, p(9 .\ ), calculations for the EM 
model (equation (9)), the MPS model (equation (11)), and the 
IMC model (equation (13)) using the LOPC data were compared 
with direct j3 measured by the 9.1 Jim CW lidar. For a given 
LOPC datum, p calculations using the three models are 
associated with the same experimental uncertainty; however, 
the difference in the p values shows the uncertainty involved 
with the way in which the aerosol composition is modeled. 
Thus the three models give a good envelope of the effect of 
changing the mixture models for estimating the complex 
refractive index and hence the scattering properties for 
inhomogeneous aerosols. Where the three models predict 
close values, the choice of the mixture model is not as sig- 
nificant; however, when the values differ, the model that best 
agrees with the lidar data would represent the best model for 
the aerosols. In addition, for a given air mass containing sev- 
eral LOPC data points, one particular model may overlap more 
with the lidar data than the other two models and thus can be 
considered as the best representation of the aerosol model for 
such an air mass. Therefore, even though the aerosol p varies 
over ~ 6 orders of magnitude and all three models track the lidar 
data fairly well, a time series comparison of LOPC modeled p 
with the lidar p gives an approximate indication of which 
model may be preferred for different air masses encountered 
during the GLOBE flights. 

LOPC data and lidar measurements, in section 5.1, are cho- 
sen from representative GLOBE flights over the Pacific Ocean 
with 6 to 8 hour periods at almost constant cruising altitude of 
about 8.0 to 8.5 km. Details of the GLOBE II flight tracks, 
with latitude and longitude, are given in Figure 1 of Cutten et 
al. [1996]. Four flights (F6, F10, FI 2, F14) shown here gave 
some very significant marine mid-tropospheric features and 
good horizontal sampling of changes in aerosol chemical and 
physical properties in large air masses over synoptic scales. 
Also, these showed some differences in the northern and 
southern hemispheric aerosol systems. In section 5.2, histo- 
grams of the aerosol p distribution are shown for both modeled 
and measured p. Section 5.3 compares LOPC data with the 10.6 
pm CW lidar data along with the ratio of /3(9.1) to j3(10.6) 
from both the lidar data sets and the LOPC modeling. 

5.1. Flight Data 

Figures for each of the flights show temporal variation of 
several parameters: (1) ambient RH data obtained by the 
DC-8 aircraft instrumentation data sets, (2) total aerosol num- 
ber density, No, integrated from LOPC Ch4 size distribution, 
(3) aerosol composition for SA, AMS, and DS in terms of the 
total integrated volumes, V$a> V A ms> anc * V ds * obtained from 
equations (l)-(5) for each sampled distribution during the 
flight, (4) modeled p(9A) results using the EM, MPS, and IMC 
models (equations (8)-(13)), and (5) calibrated CW lidar p (9.1) 
data. 

5.1.1 Hawaii to Samoa F6. This flight was deployed on 
May 20-21, 1990, from Honolulu, Hawaii, to Pago Pago, 
American Samoa, giving a crossover sampling from northern 
to southern hemisphere, encountering the ITCZ (Figure 3). 
Good qualitative agreement was obtained between the CW lidar 
p and calculated p from the three models, showing similar tem- 
poral variation. A similar trend was also seen in the micro- 



Time (UTC) 


Figure 3. Time series of F6 data from Honolulu, Hawaii, to 
Pago Pago, American Samoa, on May 20-21, 1990: (a) atmos- 
pheric relative humidity (RH), total aerosol concentration ( No , 
number of particles per cubic centimeters of air), and total vol- 
umes of each component (V54, Vams* anc * ^ ds » |ini 3 /m 3 of air) 
obtained from LOPC size distributions using equations (l)-(5) 
and integrated over all size bins; and (b) comparison of the 
p( 9.1) lidar data with the EM, MPS, and IMC models using 
equations (8)-(13). MPS model results are given as solid line 
showing the time duration of each LOPC sample, whereas EM 
and IMC are shown as diamonds and crosses, respectively, at 
the midpoint of the sample time to minimize data overlap. 


physical LOPC total No data and component volume fractions; 
however, relative volumes of each component changed during 
the flight, affecting /3. Initially, there was more V SA content, 
suggesting Hawaiian volcanic influence and possible dimethyl 
sulfide (DMS)-derived SA mixed in with a frontal development 
that was taking place near the island. This led to slightly ele- 
vated p that dropped further away from the island and 
approaching the tropical region. Around 2245 UTC the ITCZ 
was encountered with high RH, strong convective updrafts, 
variable cloud cover, and variable p. The side-viewing video 
record on board the DC-8 showed that the aircraft was flying 
mostly in a cloud-free low to moderate p region between high- 
level cirrus clouds and low-level cumulus/altocumulus clouds. 
However, occasionally, the convection would be so deep that 
the altocumulus and the cirrus would thicken and merge 
completely associated with significant increase in p. Strong 
updrafts were seen in the concurrent European Commonwealth 
Meteorological Weather Forecasting (ECMWF) data (E. W. 
McCaul Jr., private communication, 1994). The convectively 
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Figure 4. Same as Figure 3 for data from F10 from 
Christchurch, New Zealand, to Melbourne, Australia, on May 
28, 1990. 



active regions in the ITCZ lifted the high-RH marine boundary 
layer (MBL) air, possibly detraining it into the cloud-free 
regions between the cirrus and the altocumulus through cloud- 
pumping processes. Interestingly, along with enhanced total 
No , the composition of the particles (aerosol or cloud) also 
showed this effect, giving a higher fraction of nonvolatile 
component in the heated channels Ch3 and Ch5. Usually in 
mid-troposphere, Ch3 and Ch5 reflect the nonvolatile DS 
component, but in this case there was no nearby source of DS. 
Furthermore, strong upward motion registered in ECMWF data 
suggest that this large size nonvolatile aerosol content would 
most likely be nonvolatile sea-salt particles carried upward 
from the MBL through convective cloud-pumping processes 
up to the ~8 km flight altitude. Assumption of DS versus sea 
salt in the ITCZ only gives -26% variation, as the DS to sea- 
salt refractive index change is offset by the hygroscopicity 
change registered in size and high water content, lessening the 
total effect on p. Choosing sea salt in the IMC model for Ch3 
or Ch5 in the ITCZ would give P - 2.97 x 10' 8 nr 1 sr' 1 at 2307 
UTC and ft ~ 7.96 x 10' 9 m' ] sr* 1 at 2325 UTC as in comparison 
with DS (shown in Figure 3), which gave p -2.22 x 10' 8 m' 1 
sr’ 1 and 1.07 x 10' 8 m' 1 sr' 1 , respectively. Even so, the change 
in p in the ITCZ was both measured by the lidar and predicted 
by the models. Furthermore, these aerosols also contained 
more AMS component than SA. Outside the ITCZ region after 
2330 UTC, subsidence brought stable cleaner air, lowering p 
to very clean tropical air values, sometimes even below the 
detection thresholds. In the southern hemisphere a stable 
clean air mass was encountered almost to Samoa with 
comparable SA and AMS content, giving uniformly low p 
values. Lidar P and LOPC No showed a similar trend, giving 


positive correlation showing the consistency between the two 
data sets, as these quantities are directly proportional to each 
other. Also, there was occasional positive correlation between 
RH and No and between RH and p, for example, in the ITCZ, 
while at other times, for example, in tropical air masses, they 
were poorly correlated, as there is no straightforward 
relationship between the two parameters. 

In this flight, three distinct air masses were sampled: the 
subtropical air mass near Hawaii with some superimposed 
island land mass effect, an extremely convective air mass in 
the ITCZ’s with variable backscatter, and then a clean, more 
stable, tropical air mass all the way to Samoa. All three 
models agree reasonably well with the lidar p, providing a 
good envelope of modeled P for different possible scenarios of 
mixing aerosol components. Within the envelope the overall 
experimental uncertainty (LOPC size measurement errors, 
growth errors, and refractive index errors) is common to all 
three models; hence these models can be intercompared to see 
which one may be preferred for a given air mass. For tropical 
and subtropical air mass in F6 the MPS model overall was in 
closest agreement with the lidar p, but IMC and EM were not 
very far off. This finding suggests that in the tropics and 
subtropics a fairly simple aerosol system was encountered 
with possible partial mixing of SA and AMS. Since there 
wasn’t much DS in these areas, it does not matter whether or 
not one chooses a composite mixture. Agreement of the EM 
model with lidar p implies that the simple pure model works 
well for newly formed aerosols existing in clean and fresh 
conditions. Occasionally, when aerosol loading was very low 
accompanied by low p values, the LOPC modeled p was even 
lower than lidar p, with the three models giving quite different 
results. This discrepancy was caused by the LOPC’s inability 
to sample effectively the more dispersed fewer large particles 
that provide major contribution to p as compared to the lidar, 
which has larger sample volume and fewer count statistics 
problems. In addition, this sampling inefficiency in the larger 
size bins can slightly bias one channel with respect to others, 
adding uncertainty in the estimation of large-particle 
composition and v; and thereby leading to deviations in p 
results from the three models. 

5.1.2 New Zealand to Australia F10. A comparison similar 
to that for F6 was performed for F10, deployed on May 28, 
1990, from Christchurch, New Zealand, to Melbourne, 
Australia (Figure 4). This flight showed relatively low aerosol 
counts and low p ~ 10‘ 30 m' 3 sr' 1 in the southern hemispheric 
late fall clean conditions. However, there is less fine-scale 
variation than can be seen in F6 and more larger-scale vari- 
ability, suggesting that change in air masses was encountered 
during this long transit flight. There is also more variability 
in No and the component fractions. Overall, the large-particle 
concentration was low; however, the small-particle concentra- 
tion was considerable, giving moderate No values, which 
allowed reasonable differentiation of the component fractions. 
In comparison the DS component was quite low almost 
throughout this flight, showing less land mass influence in 
the southern hemisphere over the Pacific as expected. Interest- 
ingly, these small-size particles were mostly dominated by the 
SA component in southern hemisphere mid-tropospheric aero- 
sols. This domination suggests a possible stratospheric/ 
tropospheric exchange, leading to higher SA from volcanic 
origin, probably from the eruption of Kelut in February 1990, 
which dumped small-sized SA aerosols in the stratosphere, in 
agreement with the observations of Kent et al . , [1995]; how- 
ever, fresh SA may also occur from gas-to-particle conversion 




SRIVASTAVA ET AL.: AEROSOL BACKSCATTER MODELING AND LIDAR DATA 


16,613 



Figure 5. Same as Figure 3 for FI 2 data from Darwin, 
Australia, to Tokyo, Japan, on May 31, 1990. 


[Clarke, 1993]. Two cases of tropospheric/stratospheric 
exchange were seen in ECMWp data (E. W. McCaul Jr., private 
communication, 1994) just after 0300 UTC and around 0630 
UTC with lowering of the tropopause and possible 
entrainment of stratospheric air, making both No and p dip 
low. RH data were mostly anticorrelated with No, V$ A , and p. 
V SA and V A ms nearly follow each other, suggesting that the 
entire air mass sampled was of similar composition with low 
V DS , except for regions of stratospheric SA intrusion. In gen- 
eral, the southern hemispheric No was lower than the northern 
hemispheric aerosol loading but was still slightly higher than 
the clean tropical air masses encountered in F6. 

As in F6, the three aerosol mixture models show the same 
trends in p as those seen in the lidar data. Any of the three 
models could possibly be chosen for such an air mass; how- 
ever, occasionally, there is deviation between the lidar p and 
the modeled p. Even though the air was fairly clean with low p , 
the EM model (generally applicable for clean fresh aerosol 
system) seems to depart more from the lidar data than it did for 
F6, while the MPS and IMC models give more overlap. This 
finding suggests that possibly the aerosol system may not be 
as relatively fresh as in the tropics, being more aged with fine 
mode residual of either volcanic or continental origin with 
more complex mixture. 

5.1.3 Australia to Japan F12 On May 31, 1990, FI 2 
crossed over from southern to northern hemisphere, again fly- 
ing through the ITCZ in a different tropical region from 
Darwin, Australia, to Tokyo, Japan, where an Asiatic dust 
plume was encountered. Figure 5 shows the comparison for 
FI 2. Initially clean subtropical conditions followed by the 
ITCZ (around 0200 UTC) were encountered in the southern 
hemisphere. In the ITCZ, LOPC data measured broad size dis- 


tributions in all channels, suggesting large-size nonvolatile 
sea-salt aerosols yielding high p as in a similar feature 
observed in F6, associated with strong convective pumping 
carrying large-size MBL particles at the aircraft altitude. In the 
very clean tropical region near the equator (crossed around 
-0235 UTC) the aerosol loading and p dropped to threshold 
levels of both instruments. In addition, very high fine-particle 
counts were reported with a condensation nuclei (CN) counter 
on board the aircraft [Clarke, 1993], showing a region of sig- 
nificant new particle production that was mostly pure SA. Fur- 
ther on, the aerosol concentration steadily increased in the 
northern hemisphere to the levels of marine subtropical air 
mass until 0430 UTC. Even though p conditions were low, 
there was a lot of variation in comparison with the more stable 
subtropical air mass in F6. After 0430 UTC there was a distinct 
change in the air mass with higher No and more SA, AMS, and 
DS component (Figure 5a) with the encountering of the Asian 
continental landmass effect with moderate p values. Around 
0600 UTC, DS increased significantly along with p. 

Till -0430 UTC for subtropical air mass, most of the aero- 
sol models show fair agreement, though the overall best fit is 
given by the MPS model. After -0430 UTC, when the Asian 
continental aerosol was encountered with more AMS and DS, 
the IMC model gave closer overlap. This aerosol system was 
aged, probably well mixed, though diluted by clean marine 
tropospheric air until -0600 UTC. Afterward the air mass 
showed pronounced continental effect at the flight level. In 
the continental DS plume, with or without marine air dilution, 
the aerosol was most likely an internally mixed composite; 
hence the IMC model led to the closest agreement with lidar p. 

5.1.4 Japan to Hawaii F14. A transit flight FI 4 was 
deployed on June 3-4, 1990, from Tokyo, Japan, to Honolulu, 
Hawaii, giving a good sampling of northern hemispheric 
remote regions with the effect of transport of large continental 
air mass over the ocean. Figure 6 shows the data comparison. 
Initially, extremely high p due to Tokyo pollution was 
encountered superimposed on the Asian dust plume event in a 
high-RH environment. After 2400 UTC, the air mass changed 
with high DS and AMS content. Around 0200 UTC, some 
high-RH air mass was encountered for a short period 
representing cloudy air. In this case, high lidar p was 
encountered, and the modeled p showed similar increase but 
was somewhat lower, probably because of underestimation of 
D a of large wet aerosols and cloud particles in very high RH. 
After the brief high RH period, however, drier air mass than 
before was encountered. Changes in p in this air mass showed 
variations in the continental plume due to mixing and 
consequent dilution with clean marine tropospheric air. Nearer 
Hawaii at -0345 UTC the air mass abruptly changed with p 
decreasing by about 2 orders of magnitude to very clean 
subtropical mid-tropospheric levels with almost no 
continental effect. Lidar P was low, close to its detection 
threshold. The particle counts also showed a dramatic drop-off, 
with lower RH and composition changed to primarily SA with 
only occasional AMS and negligible trace DS. Finally, during 
landing in Hawaii, higher P was encountered for a short period 
due to some local AMS and DS source. 

For high aerosol loading till -0400 UTC, all three models 
gave p close to lidar data; however, the IMC model gave the 
best overlap with lidar p , again suggesting that the aerosols 
were aged well-mixed composites that have been advected over 
the clean marine mid-troposphere. For clean subtropical case 
the MPS model gave the best agreement with lidar data, while 
the IMC model, which is better suited for aged well-mixed 
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Figure 6. Same as Figure 3 for F14 data from Tokyo, Japan, to 
Honolulu, Hawaii, on June 3-4, 1990. 


aerosols, did not agree as well. This finding suggests that the 
clean subtropical aerosol system is primarily composed of 
small-size SA aerosols with partial ammoniation as shown by 
the measured LOPC V r There was some variation in the LOPC 
data, and different models showed higher uncertainty during 
very clean conditions due to the inefficiency of large-particle 
sampling as discussed earlier. During landing, all three models 
gave similar p results, well within the rapidly changing lidar 

p- 

The flights shown here sampled different mid-tropospheric 
regions and conditions during GLOBE. For some GLOBE 
flights the preferred models predicted for a given air mass have 
been used with the LOPC data and compared with the FSSP data 
calculations of p, giving good agreement [Cutten et al. 9 1996]. 

5.2. Aerosol Backscatter Distribution 

The frequency of occurrence of p{9. 1) in Figure 7 shows the 
marine mid-tropospheric aerosol backscatter distribution sam- 
pled during the four flights. Data in Figure 7a show the mode 
of p ~ 3.5 x 10" M m' 1 sr* 1 , which agrees with the mode in Fig- 
ure 6b of Cutten et al. [1996] (both of these figures are subsets 
of the overall GLOBE flight statistics shown in Figure 2 of 
Rothermel et al. [1996a]). Figure 7a has much better smooth 
statistics with about 100 times more lidar data points than the 
109 LOPC data points shown in Figures 7b, 7c, 7d for the 
three aerosol models. The LOPC data gave relatively more 
extremely low P occurrences, because in very clean atmos- 
pheric conditions it could not effectively sample the few large 
particles that would tend to increase the total p, whereas the 
lidar with its large sample volume could detect them well and 
consequently was not biased toward very low p. In spite of the 



log P(9.1 ) 

Figure 7. Comparison of the histograms of measured /3(9.1) 
from CW lidar and modeled >3(9. 1) using the LOPC data for F6, 
F10, FI 2, and F14 (shown in Figures 3, 4, 5, and 6). 

overwhelming difference in data resolution and the sampling 
volumes between the lidar and the LOPC, all three p models 
show the presence of the mid-tropospheric background aerosol 
[Rothermel et al, 1989, 1996a; Cutten et al., 1996; Pueschel 
et al, 1994] similar to the lidar p. However, the MPS model 
(Figure 7b) gave very good agreement with the lidar p, both 
showing the background mode ~3.5 x 10* 11 m' 1 sr' 1 . The IMC 
and EM models were slightly offset in the mode value, 
predicting it at — 1 x 10' 10 nr 1 sr' 1 . 

The use of the MPS model represents the first time that the 
measured aerosol microphysics has been able to predict the 
value of the aerosol background p mode in such good agree- 
ment with lidar p. This ability also reflects the possible com- 
positional mixture of the background aerosol itself 
[Srivastava et al. , 1995], suggesting that the background 
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aerosol is primarily SA with some AMS that has had sufficient 
time for the sulfates to have been mixed with each other, but 
the trace DS component is generally unmixed. The background 
mode seems due neither entirely to a very fresh aerosol system 
nor to a very aged one. If it is primarily the freshest short- 
lived pure SA, then the EM model would tend to agree; if it is a 
well-mixed, very aged coarse aerosol, then the IMC model 
would agree. Instead it seems to be a residual aerosol system 
emerging from an equilibrium between the sources (that 
produce fresh aerosols into it) and the sinks (that remove aged, 
well-mixed aerosols from it). Thus this quasi-stable 
background aerosol may predominantly consist of partially 
mixed SA and AMS, in an MPS mixture, with trace DS. 

All three models also showed another higher-/? mode at -1.7 
x 10' 9 m _1 sr _1 , comparing well with the lidar p. This type of 
analysis shows that different mixture models not only give 
different results for a single data point but also could give 
slightly different statistical results. Comparison with direct 
lidar data gives some indication of the chemical structure of 
the dominant aerosol mode in the MPS model. 

5.3, Backscatter Ratio for 9.1 and 10.6 (im 

The other CW lidar measured P at 10.6 pm, which gave a 
two-wavelength P comparison, leading to some aerosol com- 
position diagnosis [Srivastava et al., 1995]. For most of F6 
and F10, /?( 10.6) was below the lidar’s sensitivity level. 
Figure 8a shows comparison of modeled /?( 10.6) from LOPC 
data with lidar /3( 1 0.6) for F12 and F14 compiled together. 
There was very good agreement between the IMC model and 
lidar p for the aged continental plume over the Pacific, whereas 
for the fresh small-size aerosols in clean conditions the MPS 
model was in better agreement. This agreement reconfirmed 
the results of the comparison of lidar /?(9.1) and modeled 
P(9.\) (Figures 5 and 6). 

The ratio P R = P(9 . 1 )//?(10,6) is shown in Figure 8b, giving 
a comparison of the lidar Pr [Srivastava et al., 1995] with the 
p R obtained from LOPC data. Here, Pr > 2 indicates the pres- 
ence of an AMS component [Srivastava et al., 1995]. In gen- 
eral, the measured lidar p r and modeled Pr agreed very well, 
showing overall similar trends. However, for conditions of 
rapidly changing air masses (for example, -0600 UTC on the 
approach to Tokyo and -0400 UTC nearer Hawaii), LOPC data 
do not show the rapid aerosol variation detected by the lidar 
(as discussed in section 2.2). In low -p conditions 10.6 pm 
lidar data were close to the detection threshold of -2 xl O' 11 
m" 1 sr' 1 , which made p{ 10.6) more uncertain, thus adding to 
lidar Pr variability (for example, around 0300 to 0400 UTC 
near Darwin and around 0400 to 0500 UTC near Hawaii). Thus, 
in view of the major differences in the respective instrumenta- 
tion, sampling, and their limitations, comparison of the two 
methods, lidar P R and LOPC modeled P R , showed reasonably 
good agreement. 

Comparison of the inversion with lidar p R with volume frac- 
tions performed by Srivastava et al. [1995, Figure 4] with that 
of the volume fraction estimations, V' from the LOPC data is 
shown in Figure 8c. With p at two wavelengths, only two 
compositions are distinguishable. From 9.1 and 10.6 p R , AMS 
can be distinguished (because of its spectral signature at 9.1 
pm) from dominant SA in clean conditions or from DS in con- 
tinental dust, but it is harder to distinguish between SA and DS 
for intermediate cases. Hence we compared the lidar p R - inferred 
AMS volume fraction, V AMS with the more direct LOPC- 
estimated (from equations (l)-(5) and Figures 5 and 6). 


The SA and DS volume fractions combined together can be 
obtained from V$A + Vds =1 - W On average, LOPC esti- 
mates of V' AMS component showed similar trends to those 
obtained in lidar- inferred V' AMS in agreement within a factor 
of 2. 

Considering the majqr differences in the respective instru- 
mentation, sampling, limitations, and uncertainties (see sec- 
tion 2.2) that can contribute to discrepancies in the two data 
sets, agreement of Km within a factor of 2 is reasonably 
good. The p R inversion to V AMS from the remotely sensed 
lidar data is subject to the combined two -p uncertainties of 
-33% to 52% with possible uncertainties of the simple inver- 
sion method used (discussed by Srivastava et al. [1995]). 
Though the aspirated aerosol size distribution data from LOPC 
has a lower uncertainty, of -10% to 20%, it is subject to large 
particle count (major contribution to the volume) statistics 
error, which adds more uncertainty to its Vj estimation. Fur- 
ther, there appears to be a general bias in the lidar V' AMS , 
which is lower than LOPC V AMS by roughly 0.1 to 0.2. This 
is most likely due to the presence of bisulfates, which were 
included but not differentiated in LOPC V A m$ from the differ- 
ence between Ch2 and Ch3, whereas the measured lidar p R was 
inverted to V AMS , assuming purely AMS composition and 
refractive index. (Bisulfate refractive indicies are not well 
known in infrared, but preliminary absorbance data 
(S. Johnson, Argonne National Laboratory, private communi- 
cation, 1991) shows that at 9.1 pm it is lower than AMS, 
being diluted by the extra acidic component). Bisulfates would 
tend to lower the effective sulfate (AMS+bisulfate) complex 
refractive index and would theoretically predict lower /3( 9.1) 
and p R . Thus, to account for the measured p( 9.1) to /?( 10.6) 
ratio, more AMS+bisulfate component would have to be pres- 
ent (as seen in LOPC data) than if it were pure AMS (as inferred 
by lidar). So in the measured lidar Pr inversion, if an 
AMS+bisulfate mixture were to be considered, then inferred 
V' AMS would come out higher than the pure AMS case. Thus, 
under the circumstances with lidar and LOPC differences and 
the additional unaccountability of bisulfates, Pr and V’am 
comparisons do provide the best possible validation for the 
remote sensing, high-resolution, dual-wavelength lidar p R 
ratio technique. Multiple-wavelength lidar observations, set at 
other component-sensitive wavelengths, would allow further 
differentiation of SA, DS, and bisulfates. 

6. Conclusions 

Depending on the atmospheric conditions, aerosol distribu- 
tion and composition vary dramatically, affecting the scatter- 
ing properties of aerosols. Since the scattering properties of 
atmospheric aerosols, which have a major impact on the 
Earth’s climate [Charlson et al., 1991], can vary 
significantly, it is important to characterize them on a large 
scale [Penner et al., 1994]. The GLOBE flight mission 
provided large-scale sampling of aerosol microphysical, 
chemical, and backscattering properties simultaneously, under 
a wide range of atmospheric conditions. 

Component fractions inferred from the LOPC were used in 
three different mixture models: externally mixed (EM), mixed- 
phase sulfate (MPS), and internally mixed composite (IMC). 
The main conclusion of this paper is that all three models give 
fairly good agreement with the lidar p and provide a well- 
defined envelope of the modeling uncertainty itself. However, 
within this envelope under given atmospheric conditions that 
govern the complexity of the aerosols, one model sometimes 
overlaps with lidar data more than the others. The relative 
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Figure 8. Time series of (a) /J(10.6) combined for F12 and F14, showing comparison of CW lidar data with 
EM, MPS, and IMC models, (b) ratio = $9.1)/$10.6) from the lidar data in comparison with EM, MPS, 
and IMC model, and (c) the inferred volume fraction of AMS component, V^ MS from the lidar data as 
compared with the LOPC data using equations (l)-(5) and integrating over all size bins. 


choice of a model is suggested from modeled p values that are 
most consistent with the lidar fi by a comparison of several 
data points in an air mass or a distinct meteorological feature. 
The three models chosen for varying complexity of an aerosol 
system do seem to follow the complexity of the aerosols. The 
simple EM model is close to lidar data mostly for fresh undis- 
turbed aerosols. However, in general, there is more likelihood 
of partial mixing of the more reactive aerosol components 
such as SA and AMS, chemically leading to some bisulfates 
whose volume fractions, and chemical and optical properties 
are difficult to quantify. Using physically mixed SA and AMS 
in an MPS model gave the best overlap with the lidar data for 


clean mid-tropospheric aerosol systems. For more well-mixed 
aged continental-type aerosols in higher loading conditions, 
such as continental dust plumes and pollution layers, the IMC 
model gave excellent agreement with the lidar data. The 
background aerosol feature exhibited both by lidar and by 
modeled ft was best characterized by the MPS model. This also 
indicated that possible composition of background aerosol 
may be partially ammoniated SA containing trace amounts of 
DS. Thus the two lidar and LOPC data sets, keeping in view 
their individual caveats, complement each other well, giving a 
more comprehensive picture of the complex atmospheric 
aerosol system. 
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